Nodal domain distribution of rectangular drums * 
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We consider the sequence of nodal counts for eigenfunctions of the Laplace-Beltrami operator in 
two dimensional domains. It was conjectured recently that this sequence stores some information 
pertaining to the geometry of the domain, and we show explicitly that this is the case for the family 
of rectangular domains with Dirichlet boundary conditions. 
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I. INTRODUCTION 

One of the major problems in mathematical physics is 
concerned with the geometrical information stored in the 
spectrum of the Laplace Beltrami operator 



A^(r) = P^(r); r G fi(a) 



(1) 



The spectrum is ordered such that Ej-\ < Ej < Pj+i 
and f2(a) is a connected compact region, parameterized 
by a, on a 2D Riemannian manifold. If f2(a) has a 
boundary, Dirichlet boundary conditions are assumed. 
The corresponding physical system could be a vibrating 
drum. In 1911 H. Weyl showed that the number of eigen- 
values up to energy E is 



(2) 



where A is the area of Q. Subsequent research have shown 
(see e.g., pj) that each of the terms in the asymptotic se- 
ries of N(E) provides further geometrical information on 
the boundary. This prompted M. Kac to ask, 'can one 
hear the shape of a drum ?' That is, 'is it possible to 
uniquely define the shape of the drum from the spectrum 
?' It is known by now that for certain classes of domains 
the answer to Kac's question is positive, whereas there 
exists a large set of isospectral domains which are not 
isometric. (Ref. gives an updated review of this sub- 
ject.) 

In the present note we would like to investigate the ge- 
ometrical information stored in yet another sequence of 
numbers which are derived from the eigenfunctions ipj. 
Considering real eigenfunctions ipj, we count the num- 
ber Vj of nodal domains which are the connected domain 
where ipj has a constant sign. The nodal domains are 
separated by the nodal lines where ipj = 0. The sequence 
is the sequence of nodal counts. According to 
Courant's Nodal theorem Vj < j. This fundamental the- 
orem reveals the deep connection between the spectrum 
and the nodal count. It is convenient to define the nor- 
malized nodal domain numbers £j = Uj/j. Because of 
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Courant's theorem < £j < 1. This estimate has been 
further refined (for domains in R 2 ) £| 



limsup = 0.691 



(3) 



Following p| , we study the distribution of the normalized 
nodal numbers in the spectral interval I = [E°, E 1 ] 



(4) 



E 3 ei 



where Nj is the number of levels in the interval /. 

In Ref. p| the above distribution has been introduced 
as a tool to distinguish between systems which are inte- 
grable (separable) or classically chaotic. For the class of 
separable domains, it was shown that the limit distribu- 
tion 



P(0= Hm P&I) 



(5) 



exists. This has universal features: (a) there exists a 
system dependent parameter maximum value of the 
nodal domain number, such that P(£) = for £ > £' and 
(6) for £ « 



P(0 = 



C 



(6) 



The constant C is system dependent, but the order of 
the singularity is universal and depends only on the di- 
mensionality. (It was recently shown that the exponent 
for domains in d dimensions is {d — 3)/2.) 

The dependence on the geometry of the domain can 
come only through the parameters £',C or the details 
of the function P(£) away from the universal domain. 
Indeed, the limiting distributions for the rectangular and 
circular boundaries were computed in [j| and found to be 
different as expected. However, as will be shown below, 
the function P(£) does not distinguish between different 
rectangles. That is £' = 2/7T and 



-1/2 



(7) 



for all rectangles! Note that for £ w 2/ir, this result coin- 
cides with the universal expression © with C = 1/V2. 
The new result of the present note is that the depen- 
dence of P(£, I) on the finite spectral interval / contains 



2 



sufficient information to resolve between different rectan- 
gles. Thus, by counting nodal domains one can deduce 
the shape of the (rectangular) drum. It should be empha- 
sized at the outset that the nodal count sequence involves 
dimensionless integers, and therefore it cannot provide 
any scale information. Hence, when we say "resolve" we 
mean "resolve up to a scale" . 



II. RECTANGLES 

We consider the Dirichlet spectrum of a domain 
bounded in a rectangle with sides L x and L y . Denot- 
ing a = L x /Ly and choosing L x — tt, the spectrum is 
given by 



E = n 2 + a 2 m 2 



(8) 



where n,m — 1,2,3... and < a < 1. Since the sys- 
tem is separable in rectangular co-ordinates the nodal do- 
main number is simply Vj — nm, and j = N(n 2 + a 2 m 2 ) 
where N(E) is the spectral counting function. The lead- 
ing terms in the asymptotic expansion of N(E) are 



where 



Ei 2 (1 + a) 



The integral reduces to 

E° rl 



aNi J 1 



2 cos26> 



7T (1 — e/x) 



(16) 



dx (17) 



where sin 29 = %f (1 - § ) and 



I = 



if e< I 

if 2 < 2 1 



(18) 

Note that P(£,I) = for £ > £ j^. The above integral 
can be rewritten as 



2a7V/ 7 j + bx + cx 2 



where 



N(E) 



i r 

4tt 



AE- Ly/E 



(9) 



where A, L are the area and perimeter of the boundary 
respectively 0. In terms of a, 



4a I tt Je 



Introducing the transformation 
n(E,9) = V£cos(9 ; m(E,6) 



(10) 



Esm9/a, (11) 



the normalized nodal-domain number can be approxi- 
mated by 



€j(E,6) = - sin 26> 



!_£(! + «) 



(12) 



Converting the summation in eq. Q into an integral, we 
obtain the leading terms in the asymptotic expansion of 
P(£, J) in the large E limit 



1 



2aN T 



E 1 r ir/2 



E° JO 



£-£,■(£,0) dEdO (13) 



where 



^{(^-^)-l(l + a)(^-VP)}. 



Introducing the variable x = */ E/E° 



x 5 



2 sin 26* 



7T (1 — e/x) 



(14) 
(15) 



< 2 (^) , b = 2,(^ 



(f ) • (-'") 



This integral can be computed for any given value of the 
parameters 7]. 




FIG. 1: Typical behavior of the derivative P' for £ < 2/n. For 
£ > 2/tt, P is not defined as the function P is not smooth. 



III. RESULTS 

Using the above expression, it is possible to show that 
the derivative 



p , = dP 
da 



(21) 
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FIG. 2: Nodal domain distribution for the rectangular bound- 
ary. Solid, dashed and dotted curves are the approximate 
distribution (TJJ with a = 0.13, 0.44, 0.92 respectively for the 
energy range I = [10 2 ,10 4 ]. This may be compared with the 
limiting distribution \Y2'2t shown as a thick curve. 

is positive for all values of £ < 2/tt. Moreover, P', and 
hence the sensitivity to a, is maximal in the vicinity of 



the critical value f = 2/tt, as can be seen in Figure 

In Figure [3 the nodal domain distribution given by 
the eq. Ijl9(l is shown for different a, along with the 
corresponding numerical data. The limiting distribution 
is obtained by taking the spectral interval to infinity. In 
this limit, e = and 

!r 9I-V2 
L V w ; J 1 (22) 

0, i > 2/tt 

which is independent of a. Thus the parameter depen- 
dence is arising from the leading finite energy correction 
toP(CJ). 

The problem studied above shows clearly that the 
nodal sequence stores geometrical information, which, 
in the present case suffices to determine unambiguously 
the rectangle for which the nodal sequence is given. 
Attempts to generalize these ideas to other separable 
systems such as e.g., smooth surfaces of revolutions or 
flat tori are under way. 

Acknowledgments 

This research was supported in part by the Minerva 
center for complex systems and the Einstein (Minerva) 
center at the Weizmann Institute. Grants from the 
German-Israeli Foundation and the Israel Science Foun- 
dation are acknowledged. 



[1] C. Clark, "The asymptotic distribution of eigenvalues 

and eigenf unctions for elliptic boundary value problems," 

SI AM Review, vol. 9 (4), pp. 627, 1967. 
[2] M. Kac, "Can one here the shape of a drum ?," Amer. 

Math. Monthly, vol. 73, pp. 1, 1966. 
[3] S. Zeldich, "Survey of the inverse spectral problem," 

math.SP/0402356 
[4] A. Plcijel, "Remarks on Courant's nodal line theorem," 

Commun. Pure Appl. Math., vol. IX, pp. 543, 1956. 



[5] G. Blum, S. Gnutzmann and U. Smilansky, "Nodal do- 
main statistics: a criterion for quantum chaos," Phys. 
Rev. Lett, vol. 88, pp. 114101, 2002. 

[6] P.M. Morse and H. Feshbach, Methods of Theoretical 
Physics: Part I, McGraw-Hill, New York, pp. 761, 1953. 

[7] I.S. Gradshteyn and I.M. Ryzik, Table of Integrals, Series 
and Products, 6th edn., Academic Press, USA, pp. 94, 
2000. 



